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Abstract 

^ , We study implications of theories with sub-milhmeter extra dimensions and 

j-j I Mf ~ (1 — 10) TeV scale quantum gravity for neutrino physics. In these theories, 

the left-handed neutrinos as well as other standard model (SM) particles, are 
localized on a brane embedded in the bulk of large extra space. Mixing of neutrinos 
with (SM) singlet fermions propagating in the bulk is naturally suppressed by the 
volume factor Mf/Mp ~ 3 ■ 10"^^ — 3 ■ 10~^^, where Mp is the Planck mass. 
Properties of the neutrino oscillations and the resonance conversion to the bulk 
fermions are considered. We show that the resonance conversion of the electron 
neutrinos to the light bulk fermions can solve the solar neutrino problem. The 
signature of the solution is the peculiar distortion of the solar neutrino spectrum. 
The solution implies that the radius of at least one extra dimension should be in 
the range 0.06 - 0.1 mm irrespective of total number of extra dimensions. The 
corresponding modification of the Newtonian law is within the range of sensitivity 
of proposed sub-millimeter experiments, thus providing a verifiable link between 
neutrino physics and the gravity measurements. 



1 Introduction 



It has been suggested recently 0, H, H that the fundamental scale of quantum gravity, Mf, 
can be as low as few TeV. The observed weakness of gravity is the result of ( > 2) new 
space dimensions in which gravity can propagate Q. The observed (reduced) Planck scale, 
Mp = (AtcGn)^^^'^ = 3.4 • 10^^ GeV, where Gat is the Newton constant, is then related to the 
reduced Planck scale in 4 + dimensions, Mf (fundamental scale), by 

Mp = Mf^MfVN , (1) 

where Vat = L1L2....LN is the volume of the extra space, and Li is the size of the i - compact 
dimension. For definiteness we will assume that the volume has a configuration of torus in 
which case Lj = 27ri?j, where Ri {i = 1, 2, ...A^) are the radii of extra dimensions, so that 

Vn = {27: fR,R2...RN . (2) 
Using ([^) and (^) we get the constraint on the extra dimension radii: 

{2nrR,R,...R^ = . (3) 

(Notice that in some publications the factor (27r)^ is removed from this relation by redefinition 
of the fundamental scale: = {2tt)^/'^^+^^ Mf.) 

The phenomenological acceptance requires that N > 2, since for = 1 the radius would 
be of the solar system size. According to present measurement the distance above which the 
Newtonian law should not be changed is about 1 mm, and therefore 

L, = 2-KRi < 1 mm . (4) 

For N = 2 and Mf ~ TeV one gets from (|l|J^) i?i ~ -R2 ~ 0.1 mm which satisfies (^) 

Thus, in theories under consideration it is expected that the Newtonian 1/r law breaks 

down at the scales smaller than the largest extra dimension: Lmax- The experimentally most 

exiting possibility would be if L^ax ~ 1 — 10^^ mm, that is, in the range of sensitivity of 

-"^In a different context an attempt of lowering the string scale to TeV, without lowering the fundamental 
Planck scale was considered in based on an earlier observation in see also ||] for lowering the GUT 
scale. Dynamical localization of the fields on a (solitonic) brane embedded in a higher dimensional space-time 
has been suggested earlier in the field theoretic context , ||] , . For some realizations of this scenario in the 
Z?-brane context see, [||, | |lO| |. 
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proposed experiments 0]. As we will argue in this paper the same range is suggested by 
neutrino physics, namely, by a solution of the solar neutrino problem based on existence of 
new dimensions. 

Usually it is assumed that all large radii Ri are of the same order of magnitude. In such 
a case N > 2 would be well out of sensitivity of any planned sub-millimeter gravitational 
measurements. On the other hand, for N = 2 the supernova analysis pushes the lower bound 
on Mf to 30 TeV or even to 50 TeV ||12| implying that R < 0.01 mm, which is again beyond 
the planned experimental sensitivity. However, in the absence of any commonly accepted 
mechanism for stabilization of large radii^, the requirement of their equality is unjustified. 
In the present paper we will assume that radii may take arbitrary values which satisfy the 
fixed over-all volume (|]) and phenomenological constraints. In such a case the theory with 
several extra dimensions still can be subject of sub-millimeter test, while avoiding astrophysical 
and other laboratory bounds. As we will see, these bounds are sensitive to the shape of extra 
dimensions. 

In this paper we will discuss possible consequences of the high-dimensional theories for 
neutrino physics. In particular, we will suggest new high- dimensional solution of the solar 
neutrino puzzle. This solution implies that the radius of at least one extra dimension must 
be within 0.06 — 0.1 mm range. This observation relies on new high dimension mechanism of 



neutrino mass generation suggested in |Tj] which we will briefly recall. 

According to the framework elaborated in [Q, ^ all the standard model particles must 
be localized on a 3-dimensional hyper-surface ('brane') |^, H, Q embedded in the bulk of 
large extra dimensions. The same is true for any other state charged under the standard 
model group. The argument is due to the conservation of the gauge flux, which indicates 
that no state carrying a charge under a gauge field localized on the brane, can exist away 
from itp|||l||. Thus, all the particles split in two categories: those that live on the brane, and 
those which exist everywhere, 'bulk modes'. Graviton belongs to the second category. Besides 
the graviton there can be additional neutral states propagating in the bulk. In general, the 
couplings between the brane, ipbrane, and the bulk ipbuik niodes are suppressed by a volume 
factor: 

^brane'4^brane'4^bulk • (5) 



N 



^For some ideas in this direction see p^ . 



3 



According to (0) the coupling constant in (||) equals 

Mf m Mf 

and it does not depend on number of extra dimensions. It was suggested |n| to use this small 
model-independent coupling to explain the smallness of the neutrino mass. The left handed 
neutrino, ul, having weak isospin and hypercharge must reside on the brane. Thus, it can get 
a naturally small Dirac mass through the mixing with some bulk fermion and the latter can 
be interpreted as the right-handed neutrino uji: 

^H,,., . (7) 

Here H is the Higgs doublet and h is the model-dependent Yukawa coupling. After electro 
weak symmetry breaking the interaction (|^) will generate the mixing mass 

hvMf , , 



rriD 



M, 



p 



where v is the VEV of H. For Mf ~ 1 TeV and h = 1 this mass is about 5.6 ■ 10~^ eV. 

Being the bulk state, un has a whole tower of the Kaluza-Klein (KK) relatives. For extra 
dimensions they can be labeled by a set of A^ integers ni^n^-.-nisf (which determine momenta 
in extra spaces): i'ni,n2...ni^ r- Masses of these states are given by: 



nf 



Notice that the masses of the KK states are determined by the radii whereas the scale of the 
Newton law modification is given by the size of compact dimensions. 

The left handed neutrino couples with all with the same mixing mass (§). This mixing 
is possible due to the spontaneous breaking of the translational invariance in the bulk by the 
brane. 



In ref. ||T6[ a general case has been considered with possible universal Majorana mass terms 
for the bulk fermions. Neutrino masses, mixings and vacuum oscillations have been studied 
for various sizes of mass parameters. 

In this paper we continue to study the consequences of mixing of the usual neutrino with 
bulk fermions in the context elaborated in Ref. 0]. We consider both the neutrino oscilla- 
tions (in vacuum and medium) and the resonance conversion. We show that the resonance 
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conversion of the electron neutrinos to the bulk states can solve the solar neutrino problem. 
This solution implies R ~ 0.1 mm, thus giving connection between neutrino physics and sub- 
millimeter gravity measurements. We also discuss production of the KK- neutrino states in 
the Early Universe and in supernovae. 



2 Neutrino mixing with the bulk modes. UniversaUty 

Let us first assume that extra dimensions have the hierarchy of radii, so that only one extra 
dimension has radius R in sub-millimeter range and therefore only one tower of corresponding 
Kaluza-Klein modes is relevant for the low energy neutrino physics. The number and the size 
of other dimensions will be chosen to satisfy the constraint (|^). (We will comment on effects 
of two sub- millimeter dimensions in sect. 4 and 5.) 

The right handed bulk states, um, form with the left handed bulk components, Ui^, the 
Dirac mass terms which originate from the quantized internal momenta in extra dimension: 

_ n 

™n T^nRVnL + h.C, m„ = — . (10) 

n=-oo 

The mass-split is determined by 1/R. According to (|^, H) the bulk states mix with usual left 
handed neutrino (for definiteness we will consider the electron neutrino Vei) by the Dirac type 
mass terms with universal mass parameter: 

niD J^nR^eL, = ^^^^ Ri Q ■ lO^^eVh , (11) 

Mp ITeV ^ ' 

where h is the renormalized Yukawa coupling. 
The mass terms (|iy,0) can be rewritten as 



°° n 



rriDl^ORl^eL + mDj2 (^^R + ^-nR)^eL + I] p i^nR^nL " Z^-nijZ^-ni) + h.C. (12) 
n=l n=l 

Notice that z/ql decouples from the system. Introducing new states: 

1 _ 1 



and denoting by u'^^, z/^^ the orthogonal combinations we can write the mass terms in (p!0D as 



°° n 



rriDl^ORl^eL + V2mn ^ UnRl^eL + 51 (l^nRl^nL + T^nR^'nl) + ^.C. (14) 
n=l n=l ^ 
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Notice that the zero mode has smaller mixing mass with z/g than non-zero modes; the 
states z/^j^, z/^^ decouple from the rest of the system. 

Diagonalization of the mass matrix formed by the mass terms (|1^ in the limit oimoR 1 
gives (see Appendix) the mixing of neutrino with n^^ - bulk mode, Uni- 

tan6'„ ~ — — = — , (15) 
m„ n 

where 

V2hvM^ (16) 

^ Mp ^ ^ 

determines mixing with the first bulk mode. The lightest state, z/q, has the mass mo ^ nip, 
and others, z/„,: 

rrin^ — . (17) 
ri 

According to (^) the electron neutrino state can be written in terms of the mass eigenstates 

as 



^ N 

where the normalization factor equals 

N' = i+eT.-. = ^+'ie . (19) 

6 

From the phenomenological point of view the bulk modes (being the singlets of the SM 
symmetry group) can be considered as sterile neutrinos. Thus, we deal with the coupled sys- 
tem of the electron neutrino and infinite number of sterile neutrinos mixed. According to (0), 
the electron neutrino turns out to be the coherent mixture of states with increasing mass and 
decreasing mixing. 

The following comment concerning normalization is in order. The contribution to the 
normalization A^^ (pO]) from mass states with n = k,k + 1, .... (starting from number k) can 



be estimated substituting the sum by the integral: 

and for — > oo we get — > 0. Thus, due to decrease of mixing the effect of heavy states is 
suppressed. 
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In real physical processes with energy release Q only low mass part of the state (|T8]) can be 
produced. The states with n > QR do not appear. This leads to breaking of universality, that 
is, to difference of normalization of the neutrino states produced in processes with different 
Q. As follows from (^) the contribution of states with n > QR to the normalization equals 

AW) = il « 10-'»q7^ . (21) 

(For Mf = 10 TeV and R~^ = 3 • 10^'^ eV.) Since Q ^ 1/R for sub-millimeter scale the 
deviation from universality is negligible. This is not true for two extra dimensions of common 
size [0 (see later). For the same reason a change of kinematics of processes due to emission 
of states with nin ~ Q is unobservable. The probability of emission of the heavy states is 
negligible. 



3 Oscillations and Resonance Conversion 



Let us consider the vacuum oscillations of neutrino state produced as Ue (|18[) to the bulk 
modes. The state will evolve with time t as 



(t) = ^(z/o + eE;e^""^" 



1 



n=l 



n 



where the phases (/)„ equal 



ml)t 



2E 



and E is the energy of neutrinos. The survival probability of the 
equals: 



(22) 



(23) 



z/p oscillations then 



P=|(z.>e(t))P = ^ 



1 + 



n=l 



(24) 



Since 0„ oc n^, the oscillation picture consists of interference of the infinite number of modes 
with increasing frequencies oc ra^ and decreasing amplitudes oc l/n"^. For practical purpose all 
high frequency modes can be averaged, so that only a few low frequency oscillations can be 
observed depending on the energy resolution of detector. Using (^) we find the probability 
averaged over all the modes: 



1 



[1 + (7rV6)e]2 



90^ 



(25) 
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It is smaller than the two neutrino probability with the same mixing parameter ^ due to 
presence of infinite number of mixed states. In particular, for ^ ^ 1 we get P ^ 1 — (vr^/S)^^, 
whereas in the 2v case: P \ — 2^^. 

In the case when only lowest frequency mode (associated to Vi) is non- averaged, we get 
from (|2^) the survival probability 

1 



(1 + (vrV6)a' 
According to this equation the depth of oscillations equals 



V90 1^ ^ 2 



(26) 



[1 + (7rV6)e]2 ■ ^ ' 

Notice that due to presence of (practically) infinite number of the bulk modes which give 
just averaged oscillation result, the depth of oscillations of the lowest mode can not be maximal. 
Moreover, the relation between the depth and the average probability differs from the standard 
2z/ - oscillation case (see also ||T6[| ). 

Similarly one can find the probability with two non-averaged modes, etc.. 



In medium with constant density the mixing with different bulk states is modified depend- 

2 
k 



ing on m\ and the potential, V , which describes the interaction with medium 



V = Gf-^ (Ve - ) . (28) 



rriN V 2 

Here Gj? is the Fermi coupling constant, m^v is the nucleon mass, Y^. and y„ are the numbers 
of electrons and neutrons per nucleon correspondingly. The mode for which the resonance 



condition [17 



il " ^ (29) 



is fulfilled (resonance modes) will be enhanced: the effective mixing will be enhanced and the 
oscillations will proceed with large depth. The modes with lower frequencies (masses) will be 
suppressed; the high frequency modes, m\/2E ^ V , will not be modified. 

Let us consider propagation of neutrinos in medium with varying density p(r) keeping in 
mind applications to solar and supernova neutrinos. The energies of bulk states do not depend 
on density: 
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whereas for the electron neutrino we have 



He ~ V{p) . (31) 

Therefore the level crossing scheme (the {H — p) - plot which shows the dependence of the 
energies of levels He-, Hi on density) consists of infinite number of horizontal parallel lines 
(0) crossed by the electron neutrino line (0). In what follows we will concentrate on the 
case ^ <^ 1, so that m^) <^ m„ for all n, and the crossings (resonances) occur in the neutrino 
channels. The resonance density, p„, of H^ crossing with energy of n^^ bulk state: Hn = He{pn) 



equals according to (|30| , pT 



rniiriM o / N 

For small mixing (.^ <^ 1 ) the resonance layers for different bulk states (where the transitions, 
mainly, take place) are well separated: 

pn+l - Pn'> ^PnR = Pn , (33) 

n 

here ApnR is the width of the n^^- resonance layer. Therefore transformation in each resonance 
occurs independently and the interference of effects from different resonances can be neglected. 
In this case the survival probability Ue Ve after crossing of k resonances is just the product 
of the survival probabilities in each resonance: 

P = PiXP2X .... X Pfc . (34) 

Moreover, for Pj we can take the asymptotic formula which describes transition with initial 
density being much larger than the resonance density and the final density - much smaller 
than the resonance density. As the first approximation we can use the Landau- Zenner formula 



" ^ e-^- Eltl ' ^^^^ 
where EnR is the resonance energy which corresponds to maximal (initial) density Pmax in the 
region where neutrinos are produced: 

E "^""^^ . (36) 

2EGpPmaxiXe ~ \^n) 
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_ ml sin^ 2g„ p 
" 2E cos2e^ dp/ dr ^ ' 



is the adiabaticity parameter and sin^2^^„ ^ Since oc whereas sin^ 26'„, oc 

the adiabaticity parameter, /t„ = ki oc m^sin^26'„, does not depend on n for a given 
energy. Using this property we can write final expression for the survival probability as 

P ~ e-^"i^(^) , (38) 



where 



and f{E) is the step-like function 



dp/dr ' ^^^^ 



n, rJnR < -C/ < -C/n+l R 

Since high level resonances turn on at higher energies and n (x 1/E, the effect of conversion 
decreases with increase of the order of the resonance, n. Moreover, since in real situation the 
density is restricted from above and the energies of neutrinos are in certain range, only finite 
number of levels is relevant and the largest effect is due to the lowest mass resonance. 



4 Solution of the Solar Neutrino Problem 

Let us apply the results of previous section to solution of the solar neutrino problem. 

We choose the lowest (non-zero) bulk mass, mi = 1/R, in such a way that Am^ = 1 /i?^ is 
in the range of small mixing MSW solution due to conversion to sterile neutrino Ug |T9| : 

4^ = (4 - 10) • lO-'' eV^ . (41) 

This corresponds to 1/R = (2 — 3) 10^^ eV or 

R = 0.06 - 0.10 mm . (42) 

Using mass squared difference (|4l|) as well as maximal density and chemical composition of 
the sun we find from ([36| ) 

EiR ^ 0.4 0.8 MeV . (43) 
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Thus the pp-neutrinos {Epp < 0.42 MeV) do not undergo resonance conversion: Epp < Em 
(see ^O])), whereas the berylhum neutrinos {Ese = 0.86 MeV) cross the first resonance. 
(For smaller l/R"^ the pj9-neutrinos from the high energy part of their spectrum can cross the 
resonance and the flux can be partly suppressed.) 

The energies of the next resonances equal: EnR = n^Em, or explicitly: E2R = 1.6 — 3.2 
MeV, E3R = 3.6 - 7.2 MeV, E^r = 6.4 - 12.8 MeV, E5R = 10-20 MeV, Eqr = 14.4 - 28.8 
MeV. Higher resonances {n > 6) turn on at energies higher than maximal energy of the solar 
neutrino spectrum and therefore are irrelevant. The dependence of the survival probability on 
energy is shown in the fig. 1. The dips of the survival probability at the energies ~ EiR are 
due to turning the corresponding resonances. 

The effects of higher resonances lead to additional suppression of the survival probability 
in comparison with the two neutrino case. Therefore the parameter 4^^ which is equivalent to 
sin^ 29 should be smaller. We find that 

4^2 = (0.7 - 1.5) • 10"=^ (44) 

gives average suppression of the boron neutrino flux required by the SuperKamiokande results. 
According to (|), the value of ^ ( PD determines the fundamental scale: 

For small h the scale Mf can be large enough to satisfy various phenomenological bounds. 

Let us consider features of the suggested solution of the solar neutrino problem. The 
solution gives the fit of total rates in all experiments as good as usual 2z/ flavor conversion 
does: the pp-neutrino flux is unchanged or weakly suppressed, the beryllium neutrino flux can 
be strongly suppressed, whereas the boron neutrino flux is moderately suppressed and this 
latter suppression can be tuned by small variations of ^. 

Novel feature appears in distortion of the boron neutrino spectrum. As follows from fig. 1 
three resonances turn on in the energy interval accessible by SuperKamiokande {E > 5 MeV). 
The resonances lead to the wave-like modulation of the neutrino spectrum. (Sharp form (|38| , 
|0) is smeared due to integration over the production region.) The observation of such a 
regular wave structure with E oc would be an evidence of the extra dimensions. However, 
in practice this will be very difficult to realize. 

The SuperKamiokande experiment measures the energy spectrum of the recoil electrons 
from the reaction u e — u e [EH] . Integrations over the neutrino energy as well as the electron 
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Figure 1: The survival probability as the function of neutrino energy for the electron neutrino 
conversion to the bulk states in the Sun (solid line), 4^^ = 10"'^. Dot-dashed line shows the 
survival probability of the two neutrino conversion for the equivalent mixing sin^ 26 = 10~^. 
Dashed line corresponds to the survival probability of the two neutrino conversion for sin^ 29 = 
4 ■ 10^^ which gives good fit of the total rates in all experiments. 



energy of the survival probability folded with the neutrino cross-section and the electron energy 
resolution function lead to strong smearing of the distortion in the recoil spectrum. As the 
result, the electron energy spectrum will have just small positive slope (the larger the energy 
the weaker suppression) with very weak ( < 2 — 3 %) modulations. It is impossible to observe 
such a modulations with present statistics. 

Notice that relative modulations become stronger if mixing, ^, is larger than 10~^ and 
therefore suppression is stronger. This, however, requires larger original boron neutrino flux. 



The SNO experiment |2T|] will have better sensitivity to distortion of the spectrum. 

The slope of distortion of the neutrino spectrum is substantially smaller than in the case 
of conversion to one sterile neutrino (see fig. 1). In view of smearing effects due to integration 
over neutrino and electron energies (due to finite energy resolution) we can approximate the 
step-like function f{E) in ( ^0] ) by smooth function fapp{E) ^ sJeJEir. Then the smeared 
survival probability in the high energy range can be written as 



P 



(46) 



where a/E^ = n^'^p/ {R^dp/ dr^J Em). 

In the case of two large dimensions with i?i ~ i?2 ~ 0.02 — 0.03 mm (see sect. 5) the 
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number of bulk states, and therefore the number of relevant resonances increases quadratically: 
n^. (Here n ~ 5— 6 is the number of resonances in the energy range of solar neutrinos in the one 
dimension case.) Correspondingly, the approximating function fapp{E) will be proportional 
to E. As the result, ni ■ f{E) = const and the smeared survival probability will not depend 
on energy. In this case P{E) ^ const for E > Ei^ and there is no distortion of the recoil 
electron spectrum. For two different radii: R2 < Ri one can get any intermediate behaviour 
of the probability between that in (^61) and P = const. 

Common signature of both standard z/g — conversion and conversion to the bulk modes 
is the suppression of the neutral current (NC) interactions. The two can be, however, dis- 
tinguished using the following fact. In the case of the Uf. — Vg conversion there is certain 
correlation between suppression of the NC interactions and distortion of the spectrum. The 
weaker distortion the weaker suppression of the NC interactions and vice versa. In the case of 
— Hulk conversion a weak distortion can be accompanied by significant suppression of the 
NC events. This can be tested in the SNO experiment. 

No significant Day-Night asymmetry is expected due to smallness of mixing angle. 

Thus, the smeared distortion of the energy spectrum (for E > 5 MeV ) is weak or absent 
in the case of Ue — Vbuik conversion. However, in contrast to other energy independent solutions 
here pp— neutrino flux may not be suppressed, or the energy spectrum of pp- neutrinos can be 
significantly distorted. 



Notice that it is impossible to reproduce the large mixing angle MSW solution of the solar 
neutrino problem |19| in this context. Indeed, for ^ ~ 0{1) and as in (41) transitions in 
all low mass resonances are adiabatic, and therefore the survival probability has the form: 



2 



P{E) ~ P„(E) = sin^ e^^^, E^R<E< R. (47) 



n"^ 

(There are smooth transitions in the regions E ~ EnR-) For a given energy E the probability 
is determined by conversion in the nearest n^^ resonance with EnR < E. According to (|47| ) the 
probability decreases monotonously with energy, in contrast with observations. For instance, 
if P ~ 0.5 in the interval E = 0.5 ^ 2 MeV, then it will be P ~ 0.2 for E = 2 ^ 4.5 MeV, 
P ~ 0.06 for E = 4.5 ^ 8 MeV etc.. 

If ,^ > 1, for all the modes k with k"^ / {k + 1) < the resonances will be in the antineutrino 



channels and for k /{k + 1) > ^ in the neutrino channels p4 
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Note that solution of the solar neutrino problem due to the long length vacuum oscillations 
2^ into bulk modes implies Am^ ^ (0.5 — 5) ■ 10~^° eV^ and large (maximal) mixing. This 



leads to the following estimations 

4 ~ VAm2 ~ (0.7 - 2) 10-^eV, = 4 = (0-5 - 1-5) • lO^^eV (48) 

for values ^ = 0.5 — 0.7. Now the size of the extra dimension equals L = 6 — 20 cm which is 
excluded by existing tests of the Newton law. 



The approach opens however another possibility. Suppose that the radius R is small enough 
so that KK-excitations have negligible mixing with usual neutrinos. In this case the effect of 
extra dimensions is reduced to interaction with zero mode, vqr, only. Suppose that the same 
bulk field couples with two usual neutrinos: z/g and z/^ (or v^) generating the Dirac mass terms 

meD^oRVe + m^^oVaRV^,. (49) 



These terms lead to Dirac neutrino with mass rriD = + formed by uqr and the 

combination {rrieD'^e + ^^tD^'^t)/'"^!)- The orthogonal component is massless. In this way the 
z/g and Ufj, turn out to be mixed with the angle determined by tan^ = meo/^fiD- (Similar 
mechanism of mixing has been considered in Ref. ||l6l-) For Mf ~ 1 TeV and the original 



Yukawa couplings with bulk field h^, ^ ^ 1 we get rrieD ~ ^^iD ~ 10^^ eV which leads to 
Am? ~ 10~^° eV^ and maximal (large) z/g — mixing. This reproduces values of parameters 
required for the Ve ^ t^fi Just-so oscillation solution of the solar neutrino problem. The solution 
has however no generic signatures of the high dimensional theory and employs the latter as 
the source of neutrino mass only. 

5 Solar Neutrinos and Parameters of Extra Dimensions 

As follows from previous section, a solution of the solar neutrino problem via resonance con- 
version to the bulk modes implies the radius of extra dimension R = 0.06 — 0.1 mm and the 
fundamental scale Mf > 0.5 — 1 TeV. To satisfy the relation (|^) we need to introduce more 
extra dimensions. Let us assume that second large dimension exists with radius R'. From (^ 
we get 

1 {2TifM^R ,fMf\ 
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For Mf = 1 TeV both extra dimensions will have radii of the same size. In this case more bulk 
states are involved in conversion of solar neutrinos which will lead to absence of the distortion 
of the spectrum, as we have discussed in sect. 4. 

For Mf > 10 TeV the bulk states associated to 1/R' dimension do not participate in solar 
neutrino conversion, however they are relevant for other neutrino processes. Let us consider 
this in more details. Now the electron neutrino state can be written as 

= ^Ai^o + + ^ I ^ ^=i^n,k I , (51) 

where index n refers to dimension of larger radius i?, while k enumerate the bulk states from 
dimension R' . The sum over the states is divided into two parts: the first sum contains the 
states with /c = 0, that is, with small mass split only. This part corresponds to the one 
dimensional case discussed in sect. 2-4 and as we have shown in sect. 3, it does not lead to 
observable violation of universality. The second sum contains the towers of states with both 
large and small mass splits. Its contribution to the normalization of the state equals 

We can estimate the sum substituting it by the integral over n and k. Performing first 
integration over n from to oo and then over k from 1 to {QR') - the number of states which 
can be produced in the process with energy release Q, we find: 

AW)^^f5MnWfl',-(|;)'|f ln(W. (53) 

The difference of normalizations of the two states produced in processes with energy releases 
Qi and Q2 gives measure of the universality breaking: 

A21 = A(g2) - HQi) = HQ2/Qi)- (54) 

Taking = (2 — 4) ■ 10~^ as is implied by the solar neutrino data, Qi ~ 1 MeV and Q2 ~ 100 
GeV as well as R' / R < 0.1 we find A12 ~ (2 — 4) ■ 10~^ which is below present sensitivity. 



Notice, however, that for i?' ~ i? the violation can be at the level of existing bounds ||T5|. 



Let us consider the case of three extra dimensions with radii R ^ i?2,-R3- Assuming for 
simplicity the equality R2 = R3 = R', we find from (^: 



1 i2'nf/^M}JMfR 

^ . (55) 



R' Mp 
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For ~ 2.5- 10~3 eV and M/ = 10 TeV this equation gives l/R! ~ 3- 10"^ q^v. Performing 
calculations similar to those for one additional dimension we find parameter of universality 
violation 

^^^^^l^M^j (27r)3 M3 ^^^^ 

(V3 = {2n)^RR''^) which can be as low as 10~^ even for Q2 — Q1 ~ 100 GeV (in this estimation 
we used Mf = 10 TeV and h = l). 



6 Astrophysical Bounds. Atmospheric neutrinos 

The important bound on mixing of usual neutrinos with sterile neutrinos as well as with 
bulk states follows from primordial nucleosynthesis: production of new relativistic degrees of 
freedom leads to faster expansion of the Universe and to larger abundance of ^ife. There are 
two ways of production of the bulk states: (i) via oscillations and (ii) incoherently via chirality 
fiip. 

In the first case the electron neutrino produced as the coherent combination of mass eigen- 
states oscillates in medium to bulk states. Inelastic collision splits the state into active (z/g) 
and to sterile {vbuik) parts and after each collision two parts will oscillate independently. Os- 
cillations between two collisions average. Production rate is then the sum of the averaged 
oscillation effects over collisions. The condition that sterile states do not reach equilibrium 



puts stringent bound on oscillation parameters . 

The masses squared and mixing angles of bulk states (41, 0) implied by solution of the 
solar neutrino problem satisfy the following relation: 

4^-4 

Am^ . sin^ 2e„ ^ At = (4 - 8) ■ 10"^eV^ (57) 
where sin^ 2^„ = 4^^/n^. This "trajectory" in the /S.rn? — s\.v?20 - plot lies far outside the 



region excluded by primordial nucleosynthesis |^ . That is, production of the bulk neutrinos 
via oscillations is strongly suppressed. 

Let us consider incoherent production of the bulk states due to chirality flip. The produc- 
tion rate of an individual bulk neutrino (Pi) is suppressed relatively to the production rate of 
the left-handed neutrino (P^e) by the chirality - flip factor 

Pi /mz)^2 



P \ T 



(5^ 
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where the temperature T in the denominator comes from the propagator of a primarily- 
produced left handed neutrino. The multiplicity of the final bulk states is TR, so that the 
total bulk neutrino emission rate is suppressed as 

(!!!£)^fl = «(!|i) . (59) 



For the parameters implied by the solar neutrinos and T ~ 1 MeV this ratio is about 3 ■ 10~^^. 
Using (|59D we find the temperature, T^uik, at which production rate of the bulk states is 
comparable with expansion rate of the Universe: Tf,uik = ^exp ~ T'^/Mp: 



where ~ 1 MeV is the temperature of the neutrino decoupling. From ( |60| ) we find 
Tbuik ~ 200 GeV which is much above the "normalcy" temperature The modes from 
additional dimensions having smaller radii give even smaller contribution. 



The KK-neutrinos as well as the KK-gravitons produced in stars, in particular, in super- 
novae, increase the rate of star cooling [Q. No additional sources of cooling have been found 
from observations of the SN1987A which put stringent bound on production of the bulk states. 
The rate of the incoherent production of the bulk neutrinos in supernovae is suppressed by 
the same factor (pQ]). For temperature of the core of supernova T ~ 30 MeV the eq. ( pOf ) gives 
^buikf^ue ~ 3- 10~^^. Then taking into account that bulk neutrinos are emitted from the whole 
volume of the core, whereas usual neutrinos are emitted from the surface (neutrinosphere) we 
find that the luminosity in the bulk states is 5 - 6 orders of magnitude smaller than luminosity 
in neutrinos. Production of the bulk states via oscillations is strongly suppressed by matter 
effect. Matter suppression is weak or absent for bulk states with high mass: m ~ 10^ eV. 
However, their production is suppressed by very small vacuum mixing. 

So, we conclude that parameters required by solution of the solar neutrino problem satisfy 
astrophysical constraints. 



Let us comment in this connection on solution of the atmospheric neutrino problem via 
oscillations of muon neutrinos to the bulk states z/^ ^ i^buik- This solution requires smaller 
radius of the extra dimension: 1/R ^ \J Am^^„ ~ (5 — 9) ■ 10^^ eV or i? ~ (2 — 4) ■ 10^'^ mm and 
near to maximal mixing: ^ ~ 1. The latter means that the fundamental scale should be about 
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Mf ^ (10'^ TeV)/h. that is, about 3 orders of magnitude larger than that for solar neutrinos. 
(Notice that approximation ^ ^ 1 can not be used now to diagonalize the mass matrix and 
results of sect. 3 and Appendix should be corrected Still mixing of large mass bulk states 
is suppressed by factor 1/n and these states lead to finite averaged oscillation result. Only a 
few low mass states are relevant for non-averaged oscillation picture.) 

The oscillation parameters required by the solution of the atmospheric neutrino problem 
violate nucleosynthesis bound. Indeed, masses and mixing angles of the bulk modes satisfy 
now the relation: 

ml ■ sin^ ~ 4t ~ 6 ■ 10"^. (61) 



Nucleosynthesis bound reads |23 



^ 2 S-lO-^eV^ 
sm 26 

From these two equation we find that about 20 - 25 lightest bulk states are in the forbidden 
region. They turn out to be in equilibrium, whereas 1 or at most 2 are allowed. Production 
of the bulk states via oscillations can be suppressed if there is substantial (~ 10~^) leptonic 



asymmetry in the Universe 25 



Let us finally consider production of gravitons in the stars. The generic reason that saves 
the bulk gravitons from being ruled out by star cooling is the infrared-softness of the high 
dimensional gravity. On the language of four-dimensional KK modes this can be visualized 
as follows. The rate of each individual KK graviton emission in the star is suppressed by the 
universal volume factor 

(63) 



MpJ \MfJ MfVN 



The number of available final states is ~ T^Vat, so that the over-all rate is suppressed as 

According to this expression, the analysis of supernova core cooling gives for N = 2 the lower 
bound on the fundamental gravitational scale Mf about 30 TeV 0-50 TeV [|1^ . As it 
is clear from (0), the rate is determined by the value of Mj and it is insensitive to sizes 
of individual radii Ri as far as all 1/Ri < T. On the other hand, if some radii, Rk, do not 
satisfy this bound, then the corresponding KK modes can not be produced in the star and 
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the corresponding factor TRk in T^Vn of Eq. (0) has to be replaced by 1. The bottom-hne 
of this discussion is that standard constraint can be avoided if the extra dimensions have 
different radii. For instance, with one radius R ~ 0.03 mm and two others Rk > T (which can 
be reahzed even for Mf as low as several TeV's) the rate becomes 



/ T 



2 



r,™. ~ [j^j (TR) ^ . (65) 

Therefore, some of the radii can be of sub-millimeter size and, thus, can be a subject of direct 



experimental search in future gravitational measurements [11 



7 Conclusions 

We have studied consequences of the neutrino mixing with fermions propagating in the bulk in 
the context of theories with large extra dimensions. The bulk fermions could be components 
of bulk gravitino or other singlets of the SM gauge group. 

Phenomenology of this mixing is determined by the following features: (i) The bulk 
fermions can be considered as sterile neutrinos, (ii) Large number of these sterile neutri- 
nos is involved in physical processes, (iii) For mo < 1/R usual neutrinos are combinations of 
mass eigenstates with increasing masses and decreasing admixtures. 

The effect of bulk states with large masses is reduced to averaged oscillations. Low modes 
can lead to non-averaged oscillations in vacuum (uniform medium) or to multi-resonance 
conversion in medium with varying density. 

The resonance conversion of the electron neutrino to the bulk states can solve the solar 
neutrino problem. Properties of this solution are similar to those due to conversion to sterile 
neutrino. The important difference is that significant suppression of the boron neutrino flux 
can be accompanied by weak distortion of the energy spectrum. Moreover, weak modulation 
of the boron neutrino spectrum is expected due to conversion to several KK-states. 

Simultaneous explanation of the atmospheric, solar and LSND results in terms of neutrino 
oscillation/conversion implies existence of sterile neutrino. Moreover, the data favour ^ Vr 
oscillations as the solution of the atmospheric neutrino problem, so that the solar neutrinos 
should be converted to sterile states. In this connection one can consider the following pos- 
sibility. There is some usual (4 dimensional) mechanism of generation of the active neutrino 
masses. This mechanism produces neutrino mass pattern with heavy and strongly mixed 
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and Vr and very light z/g. Such a pattern can explain the atmospheric neutrino and LSND 
results. Neutrinos (in general of all flavors) couple also with the bulk fermion. These couplings 
(being of the same order for all neutrino species) generate negligible mixing of the KK-states 
with heavy and z/^ and large enough mixing with light Ve-, so that the solar neutrino problem 
can be solved as is described in sect. 4. 



The suggested solution of the solar neutrino problem implies that the radius of at least one 
extra dimension is in the range 0.06 - 0.10 mm, that is, in the range of sensitivity of proposed 
gravitational measurement. The fundamental scale should be about hMf ~ 1 TeV. This mass 
satisfies the fixed overall volume condition provided additional large extra dimensions exist. 
In the case of one additional extra dimension its radius should be 1/-R' = 1.3- 10~^(Mj/lTeV)^ 
eV. For h 1 one has Mj ~ 1 TeV, and ~ 2-10^'^ eV, so that the second extra dimension 
will influence the solar neutrino data too. If /i ^ 1, the fundamental scale can be much larger 
than 1 TeV, and R' can be much smaller than 1 mm. For h = 0.1 and Mf = 10 TeV we get 
1/i?' ~ 10 eV. For two additional dimensions the common radius equals 1/i?' ~ 3 ■ 10^^ GeV, 
if /i = 1 and Mf = 10 TeV. 

For large fundamental scales {Mf > 10 — 20 TeV), direct laboratory searches at high 
energies will be practically impossible and neutrinos can give unique (complementary to grav- 
itational measurements) opportunity to probe the effects of large extra dimensions. 



Appendix 



The mass terms (|1^) can be written as 



vlMur , 

where ^1 = {ul, i>iL, t>2L---) and z/^ = (z/qr, t'lij, i>2R---)] the modes z/ql, ^'ul Kr decouple from 
the system, and the mass matrix M for k + 1 states equals 



M 



( rrijj -\/2^D \f2mj:, 



V 



R 




2_ 

R 








R I 
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Let us consider the matrix MM ^ which determines mixing of the left handed neutrinos: 



MM'f 



1 



( (k + i/2)e 




2^ • 


• k^ \ 




1 


. 


. 







4 . 


. 







. 





(Al) 



where ^ = \^mj:,R. Notice that the mass matrix (Al) corresponds to compactification on the 
circle; it can be shown that the same matrix follows from the Z2 - orbifold compactification 
16| . The diagonalization of matrix can be performed starting from the heaviest state k ^ 1. 
Then one can check that the limit k ^ 00 does not change results. The rotation by the angle 
9k in the plane i/ql ~ ^feL 



4- 



{A2) 



diagonalizes corresponding submatrix. It is easy to perform the diagonalization using ^ as an 
expansion parameter: ^ ^ 1 as is implied by the solution of the solar neutrino problem. The 
rotation (Al) leads to modification of the first diagonal element {k + 1/2)^^ —>■ {k — 1/2).^^ + 
0{C,^) and modification of the mixing terms as cosdknC,. For small C,'- tanOk ~ C,/k 

and the eigenvalues equal ml ^ k'^/R'^. After k — 1 subsequent rotations we get for the first 
diagonal element: 3/2^^ + 0{^'^) and for off-diagonal term: 



^ cos 6^2 cos 6'3 . . . cos 6'/f 



« - 1^' 



E 1/ 



n 



n=2...k 



These results can be obtained from the exact characteristic equation for the mass eigenstates: 
Det [MM"^ — m^] which can be written explicitly as 

TT^^ cot(7rmi?) = 2mR. 

which coincides with the characteristic equation in |T^ for the same neutrino system. 
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